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Proper  water  management  in  a  hydrogen-fueled  polymer  electrolyte  membrane  (PEM)  fuel  cell  is  critical 
for  performance  and  durability.  A  mathematical  model  has  been  developed  to  elucidate  the  effect  of 
thermal  conductivity  and  water  vapor  diffusion  coefficient  in  the  gas  diffusion  layers  (GDLs).  The  fraction 
of  product  water  removed  in  the  vapor  phase  through  the  GDL  as  a  function  of  GDL  properties/set  of 
material  and  component  parameters  and  operating  conditions  has  been  calculated.  The  current  model 
enables  identification  of  conditions  wherein  condensation  occurs  in  each  GDL  component.  The  model 
predicts  the  temperature  gradient  across  various  components  of  a  PEM  fuel  cell,  providing  insight  into 
the  overall  mechanism  of  water  transport  in  a  given  cell  design.  The  water  condensation  conditions  and 
transport  mode  in  the  GDL  components  depend  on  the  combination  of  water  vapor  diffusion  coefficients 
and  thermal  conductivities  of  the  GDL  components.  Different  types  of  GDLs  and  water  transport  scenarios 
are  defined  in  this  work,  based  on  water  condensation  in  the  GDL  and  fraction  of  water  that  the  GDL 
removes  through  the  vapor  phase,  respectively. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Optimal  water  management  is  important  for  polymer  electrolyte 
membrane  (PEM)  fuel  cell  performance,  durability  and  rapid  start¬ 
up  under  frozen  conditions.  A  low  rate  of  water  removal  from  the 
cathode  impedes  oxygen  transport;  while  excessive  water  removal 
from  the  electrode  may  cause  electrode  and  membrane  dry-out 
resulting  in  low  proton  conductivity  and  low  apparent  exchange 
current,  as  well  as  high  losses  associated  with  membrane  resis¬ 
tance. 

The  gas  diffusion  layer  (GDL)  is  a  crucial  element  that  impacts 
the  overall  water  transport  in  a  PEM  fuel  cell.  Fundamental  under¬ 
standing  of  liquid  water  transport  in  the  GDL  and  associated 
strategies  to  keep  the  bulk  volume  available  for  oxygen  diffusion 
is  necessary.  Experimental  determination  of  capillary  pressure  as 
a  function  of  water  saturation  combined  with  Darcy’s  law  have 
been  the  primary  focal  points  in  water  management  [1-6].  How¬ 
ever,  recent  experimental  and  modeling  results  indicate  that  the 
transport  of  liquid  water  in  GDLs  is  a  process  of  capillary  finger¬ 
ing  [7-11  ],  which  cannot  be  adequately  represented  by  Darcy’s  law 
according  to  [12].  Liquid  water  flows  through  GDL  in  the  form  of 
connected  clusters,  encountering  dead  ends  due  to  the  presence  of 
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varied  diameter  pores,  and  eventually  percolates  through  a  path¬ 
way  of  the  least  resistance.  On  the  other  hand,  Owejan  et  al.  [13] 
recently  found  that  vapor  diffusion  is  an  important  mechanism  of 
water  removal  from  the  cathode  catalyst  layer  through  the  GDLs. 
It  is  shown  that  the  thermal  gradient  in  the  cathode  GDL  can  be 
sufficient  to  remove  product  water  in  vapor  state  through  the  GDL 
at  high  current  densities  [13]. 

Thermal  properties  of  the  GDL  and  its  design  are  thus  very 
important  for  high  performance  and  durability  [14-17].  Normally, 
in  an  operating  fuel  cell  there  is  3-5  °C  temperature  gradient  at 
lAcm-2  between  the  catalyst  layers  and  gas  channels  [18-20]. 
That  corresponds  to  8-19%  A RH  at  75  °C  cell  temperature.  The 
dissipation  of  the  cathode  reaction  heat  relies  on  GDL  effective 
thermal  conductivity,  which  strongly  depends  on  the  anisotropic 
packing  and  thermal  contact  of  GDL  carbon  particles  and  fibers. 
Interaction  between  mass  and  thermal  transport  makes  favor¬ 
able  water  management  a  challenge.  Only  a  handful  of  studies 
have  attempted  to  investigate  the  thermal  impact  of  GDL  on  mass 
transport.  Pasaogullari  et  al.  [14]  studied  the  anisotropic  heat 
and  water  transport  in  a  cathode  GDL.  Weber  and  Newman  [20] 
investigated  the  non-isothermal  phenomena  and  water  transport 
through  the  membrane.  They  showed  that  the  temperature  gra¬ 
dient  within  a  fuel  cell  is  sufficient  to  cause  a  heat-pipe  effect 
where  the  water  is  moving  from  cathode  to  anode  down  the 
temperature  gradient.  However,  these  studies  did  not  include  a 
micro-porous  layer  (MiPL),  which  is  a  key  component  to  the  GDL 
(see  Fig.  1). 
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Nomenclature 

D  water  vapor  diffusion  coefficient  in  porous  medium 

(cm2  s-1) 

Def f  effective  water  vapor  diffusion  coefficient  in  sub¬ 

strate  +  MiPL  composite  (cm2  s-1 ) 

F  Faraday’s  constant,  96487  (C  mol-1 ) 

AH  water  latent  heat  (J  g-1 ) 

i  current  density  (A  cm-2 ) 

L  thickness  (cm) 

P0  saturated  water  vapor  pressure  at  the  WPT  coolant 

temperature  (din cm-2) 

Qw_>v  latent  water  heat  per  unit  current  density 

(W(Acm)-2) 

specific  reaction  heat  per  unit  current  density 
(WA-1) 

R  ideal  gas  constant,  8.7  x  107  erg  (mol  K)-1 

RH  relative  humidity 

T  temperature  (I<) 

T0  WTP  coolant  temperature  (K) 

5  specific  reaction  entropy  (J  (mol  K)-1 ) 

a  fraction  of  heat  that  is  removed  by  the  cathode-side 

WTP 

r]  overpotential  (V) 

X  thermal  conductivity  (W (cm  K)-1 ) 

Aeff  effective  thermal  conductivity  of  substrate  +  MiPL 

composite  (W  (cm  I<)-1 ) 

§  fraction  of  water  removed  through  the  vapor  phase 

Subscripts 
A  anode 

SM  interface  between  the  substrate  and  MiPL 

C  cathode 

S  substrate 

M  membrane 


Although  isothermal  modeling  of  MiPLs  are  available  in  the 
literature  [21-23],  the  conclusions  from  different  groups  are  con¬ 
tradictory.  Weber  and  Newman’s  modeling  results  show  that  the 
MiPL  increases  the  back-diffusion  rate  of  water  from  the  cathode 
through  the  membrane  to  the  anode  [22];  however,  Passaogullari 
et  al.  [21]  and  Nam  andKaviany  [23]  showed  that  the  MiPL  enhances 
water  removal  rate  to  the  cathode  gas  channel.  Recent  experimental 
results  from  Karan  et  al.  [24]  show  no  statistical  impact  of  MiPL  on 
net  water  transport  across  the  membrane.  Clearly,  the  role  of  a  MiPL 
on  water  transport  remains  unresolved  (see  [24]  and  references 
therein). 

In  this  work,  thermal  effects  are  included  in  a  model  based 
on  analysis  of  a  composite  GDL  (with  a  MiPL)  to  elucidate  the 
water  transport  mechanism  in  a  porous-plate  PEM  fuel  cell.  The 
variations  in  GDL  properties  considered  are  thermal  conductivity 
and  water  vapor  diffusion  coefficients.  Experimentally  determined 
thermal  conductivities  are  used  for  quantitative  prediction  of  the 
transport  fluxes.  This  work  also  aims  to  provide  insight  into  the 
role  of  MiPL  in  water  transport  and  a  new  physical  model  to 
clarify  the  importance  of  an  intermediate  porous  layer  to  aid 
water  vapor  diffusive  transport,  as  discussed  by  Owejan  et  al. 
[13]. 

This  model  accounts  for  product  water  vapor  transport  across 
the  MiPL  present  in  a  cathode  composite  GDL  used  for  either  a  con¬ 
ventional  solid  plate  or  a  porous  plate  system  like  that  described 
in  this  work.  Porous  bipolar  plates  (as  shown  in  Fig.  1),  known 
as  water-transport  plates  (WTPs)  [25-29],  were  developed  by  UTC 
Power.  The  model  does  not  describe  liquid  water  removal  from  the 


substrate  layer  of  the  composite  GDL,  which  will  be  addressed  in 
our  future  work. 

2.  Physical  model 

2.1.  Physics  assumption 

The  WTPs  perform  two  main  functions  in  terms  of  water  man¬ 
agement.  First,  when  there  is  excess  water,  the  WTPs  provide  a 
removal  path  for  liquid  water  to  prevent  flooding.  Second,  when 
the  gas  streams  are  not  saturated,  the  WTPs  provide  water  for 
evaporation  into  the  gas  channels  and  humidifies  the  membrane 
electrode  assembly  [27].  The  liquid  water  in  the  GDL  is  comprised 
of  water  streams  flowing  from  the  cathode  to  WTP  and  pendant  or 
“stranded”  water  (as  shown  in  Fig.  1)  not  connected  to  the  water 
streams.  In  order  to  maximize  oxygen  transport,  it  is  desirable  to 
remove  the  major  fraction  of  water  in  the  vapor  phase  to  minimize 
the  fraction  of  the  GDL  volume  occupied  by  liquid  water  streams 
and  to  prevent  vapor  condensation  and  accumulation  of  pendant 
water  in  the  GDL. 

The  driving  force  for  vapor  transport  through  the  GDL  is  the  gra¬ 
dient  of  water  vapor  partial  pressure  ( WVPP).  The  gradient  of  WVPP 
is  caused  by  a  small  temperature  difference  between  the  cathode 
and  WTP,  resulting  in  the  difference  between  the  equilibrium  WVPP 
in  the  WTP  gas  channel  and  that  in  the  cathode  gas  pores.  Model¬ 
ing  results  and  experimental  data  indicate  that  in  the  present  WTP 
design,  dry  inlet  gas  becomes  fully  saturated  by  liquid  water  com¬ 
ing  through  the  WTP  at  the  gas  channel  inlet  and  stays  humidified 
along  the  entire  channel  length.  For  the  boundary  conditions  of  this 
model,  we  assume  that  water  vapor  is  in  equilibrium  with  liquid  at 
the  cathode/MiPL  interface  (x  =  1  in  Fig.  1)  and  the  substrate /WTP 
interface  (x  =  0  in  Fig.  1 ).  One  would  expect  that  there  is  zero  flux 
of  water  from  cathode  to  WTP  because  at  both  locations  water 
vapor  is  equilibrated  with  liquid  water.  Flowever,  the  water  vapor 
concentration  is  higher  at  the  cathode/MiPL  interface  than  that  in 
WTP  because  the  temperature  of  the  cathode  electrode  is  higher 
than  that  of  the  WTP.  That  concentration  gradient  creates  the  dif¬ 
fusion  flux  (see  Appendix  A  for  more  details  on  thermodynamics 
of  the  transport  process).  We  also  assume  that  no  condensation  or 


Fig.  1.  Schematic  of  a  half  fuel  cell  with  water  transport  plate  (WTP).  The  arrows  indi¬ 
cate  the  flow  of  liquid  water  in  the  WTP.  S,  M,  and  SM  represent  substrate,  MiPL,  and 
the  interface  between  substrate  and  MiPL.  Water  vapor  concentration  and  tempera¬ 
ture  at  WTP/substrate  interface,  and,  and  at  MiPL/cathode  interface,  and,  are  shown 
in  the  figure.  Selected  area  in  the  substrate  is  enlarged  to  show  the  continuous  flow 
streams  and  pendant  water. 


Table  1 

Possible  Scenarios  of  water  removal  through  the  GDL. 
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GDL  type 

Scenario  of  water 

Condensation  in  the  GDL 

Fraction  of  water  removed 

Liquid  water  streams 

Drops  and  puddles 

dynamics 

through  the  gas  phase 

in  the  GDL 

in  the  GDL 

Type  1 

Scenario  1(a) 

No  condensation  under-saturated 

100% 

No 

No 

Scenario  2(a) 

vapor  in  the  GDL 

Less  than  100% 

Yesa 

No 

Scenario  1(a) 

100% 

No 

No 

Type  2 

Scenario  1(b) 

Condensation  in  the  substrate 
Under-saturated  vapor  in  the  MiPL 

Less  than  100% 

No 

Yes  in  the  substrate 

Scenario  2(b) 

Condensation  in  the  substrate 

Less  than  100% 

Yes 

Yes  in  the  substrate 

Saturated  vapor  in  the  MiPL 

a  Note:  The  water  stream  is  coming  from  the  MiPL/cathode  interface.  There  is  only  evaporation  in  the  GDL 


evaporation  occurs  within  the  GDL;  therefore,  there  is  no  source  or 
sink  term  for  the  thermal  and  mass  transport  equation. 

The  GDL  type  defines  if  condensation  and  consequent  flooding 
occur  in  one  or  more  GDL  components.  The  scenario  of  water  trans¬ 
port  (see  Table  1)  defined  in  this  work  determines  the  fraction  of 
water  removed  through  the  gas  phase.  It  depends  on  the  combina¬ 
tion  of  water  vapor  diffusion  coefficients,  thermal  conductivities  of 
GDL  components  and  also  on  the  boundary  condition  at  x  =  0  and 
x  =  1  (see  Fig.  1). 


c=  Ls 

Ds  +  AsAm)Dm 

0) 

^  _  (AAm)Dm 
f'S  +  AsAm)£m 

(10) 

Note  that  a  +  b  =  1  and  c  +  d  =  1.  Here  L  is  the  thickness  of  each  layer 
and  the  subscripts  0,  S,  M,  C,  and  SM  refer  to  WTP,  substrate,  MiPL, 
cathode  catalyst  layer  and  substrate/MiPL  interface,  respectively 
(see  Fig.  1). 


2.2.  Mathematical  formulation 

2.2 A.  Vapor  transport  through  cathode  GDL 

In  this  section,  the  governing  equations  of  water  vapor  and  ther¬ 
mal  transport  are  introduced.  Any  in-plan  temperature  differential 
is  neglected  because  the  in-plan  thermal  conductivity  is  one  order 
of  magnitude  higher  than  that  in  through-plan  direction  [14].  Due 
to  the  uniformity  of  water  vapor  transport  in  the  PEM  fuel  cell 
plan-form,  the  transport  equations  can  be  represented  in  1-D,  in 
the  through  plane  of  the  GDL.  The  mass  and  thermal  transport 
governing  equations  are 

»*(£)-» 

*£(£)-• 

where  x  is  the  through  plane  direction  of  the  GDL,  D  the  water 
vapor  diffusion  coefficient,  C  the  water  vapor  concentration,  and 
A  is  the  thermal  conductivity.  Gas  diffusion  coefficients  in  the  sub¬ 
strates  and  MiPLs  are  corrected  by  porosity  using  the  Bruggman 
correlation.  Both  D  and  A  are  assumed  to  be  constant  in  each  layer. 

At  the  substrate  and  MiPL  interface,  fluxes  of  mass  and  heat  are 
continuous.  This  is  reflected  by  mass  and  heat  conservation  equa- 


tions  (3)  and  (4) 

CSm  -  Co  Cc  -  CSm 

d  s  T  =  D  Mr 

Ls  Lm 

(3) 

i  DSm  -  D0  i  Dc  -  Tsm 

As  T  =  am  t 

Ls  Lm 

(4) 

Solving  Eqs.  (3)  and  (4)  for  vapor  concentration  at  substrate/MiPL, 
CSm,  and  for  the  temperature  at  substrate/MiPL  interface,  Ts m,  we 
obtain 


Csm  =  aCc  +  bC0 

(5) 

Tsm  =  cTc  +  dT0 

(6) 

where 

Ds 

Ds  +  (Ds/Dm)Lm 

(7) 

k  _  (Ds/Dm)Lm 

Ds  +  (Ds/Dm)Lm 

(8) 

2.2.2.  Types  of  GDL  and  water  transport  mode 

In  this  section  we  define  different  types  of  GDL  as  shown  in 
Table  1.  The  condensation  in  a  PEM  fuel  cell  component  occurs 
when  the  local  water  vapor  concentration  at  location  x,  C(x),  is 
higher  than  the  equilibrium  concentration  at  local  temperature, 
T(x).  For  common  GDL  materials,  the  estimated  temperature  differ¬ 
ence  between  two  GDL  interfaces  is  1-3  °C  at  1 A  cm-2  (calculated 
below).  Within  such  a  small  temperature  range,  a  linear  relation 
between  the  equilibrium  vapor  concentration  and  temperature  can 
be  reasonably  assumed 

Ceq(T)  =  A  +  BT  (11) 

It  is  assumed  that  vapor  concentration  at  MiPL/cathode  inter¬ 
face,  Cc,  and  vapor  concentration  at  substrate/WTP  interface,  Co, 
are  in  equilibrium  with  liquid  water  according  to  Eq.  (11 ) 

C0(T)  =  A  +  BT0  (12) 

Cc(T)=A  +  BTc  (13) 

Substituting  C0  and  Cc  in  Eq.  (5)  with  Eqs.  (12)  and  (13)  gives 

CSM  =A  +  B(aTc  +  bT0)  (14) 

The  condition  a  =  c,  which  equals  to  b  =  d,  separates  different  GDL 
types.  That  means  water  vapor  is  either  condensed  or  under¬ 
saturated  at  two  side  of  the  condition.  If  a  =  c  or  b  =  d,  combining 
with  Eq.  (6),  the  above  equation  becomes 


Csm  =A  +  BTSm 


(15) 


According  to  Eq.  (11 ),  CSm,  defined  in  Eq.  (15),  equals  the  equilibrium 
water  vapor  concentration.  Therefore,  water  vapor  is  in  equilibrium 
with  liquid  water  at  the  substrate/MiPL  interface. 

To  maintain  a  =  c  or  b  =  d,  we  must  have 

If  the  water  vapor  diffusion  rate  is  larger  than  the  thermal  conduc¬ 
tivity,  i.e., 


Ds  >  Xs_ 
Dm  Am 


(17) 


then 

Csm  is  lower  than  the  equilibrium  concentration,  C|JJj,  deter¬ 
mined  by  Eq.  (11). 
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wvpp  (a)  wvpp  (b) 


- Calculated  vapor  pressure  in  GDI: - Saturated  vapor  pressure 


- Vapor  pressure  in  GDL  - Saturated  vapor  pressure 


Fig.  2.  The  WVPP  distribution  (dashed  red  line)  calculated  from  the  vapor  diffusion 
equations  is  compared  to  the  WVPP  distribution  of  the  saturated  vapor  (solid  blue 
line).  The  GDL  belongs  to  Type  1  if  the  red  line  is  located  below  the  blue  line  (a) 
and  it  belongs  to  Type  2  otherwise  (b).  The  diffusion  equations  are  solved  under  the 
condition  that  temperature  of  the  cathode/GDL  (x  =  1)  interface  is  higher  than  that 
of  the  GDL/WTP  interface  ( x  =  0)  and  the  vapor  is  saturated  at  these  boundaries.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  the  article.) 


Thus  in  this  case  the  water  vapor  is  under-saturated  at  the  inter¬ 
face  between  the  MiPL  and  substrate.  Two  types  of  GDLs  are  defined 
based  on  the  water  saturation  level  at  the  Mi  PL/ substrate  interface. 
The  GDL  belongs  to  Type  1  if  Eq.  (17)  holds.  Otherwise,  the  GDL 
belongs  to  Type  2.  As  shown  in  Fig.  2(a),  the  GDL  belongs  to  Type 
1  if  the  vapor  diffusion  equation  predicts  under-saturated  vapor 
in  the  GDL.  Otherwise,  the  GDL  belongs  to  Type  2,  as  shown  in 
Fig.  2(b). 

In  the  above  formulation,  we  assume  a  saturated  water  vapor 
concentration  at  the  MiPL/electrode  interface.  If  the  water  vapor  is 
under-saturated  at  this  interface,  the  water  vapor  concentration  in 
the  GDL  bulk  is  even  lower  than  that  calculated  above.  Thus,  Eq. 
(17)  is  a  conservative  criterion  for  not  having  condensation  in  the 
GDL.  All  the  GDLs  within  Type  1  have  under-saturation  water  vapor 
in  GDL  bulk.  Type  2  GDLs  would  have  either  saturated  or  under¬ 
saturated  water  vapor  in  the  GDL  bulk,  as  shown  in  Table  1  and 
Fig.  4. 


2.3.  Water  transport  scenarios 

We  are  not  only  interested  in  the  water  saturation  in  the  GDL,  but 
also  the  water  vapor  removal  rate  in  an  operating  fuel  cell.  Below 
we  analyze  water  removal  modes  through  the  GDL  in  an  actual  fuel 
cell  and  classify  them  into  different  scenarios.  The  basis  for  different 
scenarios  is  the  fraction  of  water  that  the  GDL  removes  through  the 
vapor  phase.  To  calculate  this  fraction,  we  compare  the  water  gen¬ 
eration  rate  to  the  maximum  vapor  flux  through  the  GDL  given  that 
the  water  vapor  is  in  equilibrium  with  liquid  water  at  the  GDL/WTP 
interface  and  that  the  temperature  difference  between  the  WTP  and 
the  cathode  is  fixed. 

From  Eq.  (3),  the  vapor  flux  is 


Cc-Cp 

(Ls/Ds)  +  (Lm/Dm) 


(18) 


Fig.  3.  The  WVPP  distribution  in  the  GDL  of  Type  1  is  compared  to  the  saturated 
WVPP  distribution.  The  vapor  in  the  GDL  is  under-saturated  and  there  is  no  con¬ 
densation  in  any  GDL  component.  In  Scenario  1(a),  the  vapor  at  the  cathode/GDL 
interface  is  under-saturated  and  100%  of  water  generated  in  the  cathode  can  be 
removed  through  the  vapor  phase.  In  Scenario  2(a),  the  vapor  at  the  cathode/GDL 
interface  is  saturated  and  less  than  100%  of  water  generated  in  the  cathode  can  be 
removed  through  the  vapor  phase. 


Using  the  ideal  gas  law  to  represent  the  concentration,  Eq.  (18) 
can  be  written  as 


(Pc/Tc)-(Po/T0) 
R  (Ls/Ds)  +  (Lm/Dm) 


where  P  is  the  water  vapor  partial  pressure.  Similar  to  Section  2.2.2, 
a  linear  relation  between  the  water  vapor  pressure  and  temperature 
is  assumed 


Pc  =  Po  +  (Pc  -  Po 

Then  combining  Eqs.  (19)  and  (20),  eliminating  Pc  gives 

=  (Pc-Po)(dPo/dP)-(Po/Po) 

V  PPc  (VP>s)  +  (Pm/P>m) 

Similar  to  Eq.  (19),  the  heat  flux  is 
Pc-Po 

H  (PsAs)  +  (PmAm) 


(20) 


(21) 


(22) 


Expressing  (Tc  -  Po)  and  substituting  Eq.  (21 )  into  Eq.  (22)  gives 


Qh((^s/D  s)  +  Am  Am))  ( dPp/dT )  -  (P0/T0) 
PPc  (Ls/Ds)  +  (Lm/Dm) 


(23) 


Assume  the  heat  flux  through  the  cathode  GDL  is  proportional 
to  the  total  heat  generated  in  a  fuel  cell  at  specific  current  density, 
i,  i.e., 


9h(Q  —  <^(Qr  —  Qyy^V  £)*' 


(24) 


where  §  is  the  fraction  of  water  removed  in  the  vapor  phase  and  a 
is  the  fraction  of  heat  that  is  removed  from  the  cathode  side  (our 
calculations  show  that  typically,  a  «  0.25-0.5).  The  specific  reaction 
heat  Q£  and  the  latent  heat  of  vaporization  Q°^v  are  calculated  as 
functions  of  the  specific  reaction  entropy  AS  and  overpotential  p  in 
the  following  equations 


Qg  =  T^+T^  +  q 


Qw-v  = 


2  F  4  F 

AHx  18 


2F 


(25) 

(26) 


Table  2 

Criteria  determining  the  scenario  of  water  removal  through  the  substrate. 


GDL  type 

Type  1 

Type  2 

Scenario 

Criteria  (17) 

Criteria  (27)  for  the  substrate 

Criteria  (27)  for  the  MiPL 

Criteria  (27)  for  the  whole  GDL 

Scenario  1(a) 

True 

N/A 

N/A 

True 

Scenario  2(a) 

True 

N/A 

N/A 

False 

Scenario  1(a) 

False 

True 

True 

N/A 

Scenario  1(b) 

False 

False 

True 

N/A 

Scenario  2(b) 
False 

False 

False 

N/A 
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WVPP  VWPP  VWPP 


- Vapor  pressure  in  GDL  - Saturated  vapor  pressure 


Fig.  4.  Actual  WVPP  distribution  is  compared  to  the  saturation  vapor  pressure  in  the  Type  2  GDL  In  Scenario  1(a),  the  actual  WVPP  is  lower  than  the  saturated  vapor  pressure. 
In  Scenario  1(b),  the  actual  WVPP  is  equal  to  the  saturated  vapor  pressure  in  the  substrate.  In  Scenario  2(b)  the  actual  WVPP  is  equal  to  the  saturated  vapor  pressure  in  the 
substrate  and  MiPL. 


Here  AH  is  the  latent  heat  of  vaporization  for  water.  The  anode 
and  cathode  entropy  change,  A Sa  and  A Sc,  are  130J(molI<)_1  and 
65J(molK)_1,  respectively  [30].  At  1 A  cm-2  and  0.7  V  cell  voltage, 
Q°,^v  =  0.21  WA-1  andQ°  =  0.77  WA'1.  The  temperature  profile 
within  the  cathode  can  thus  be  calculated  using  available  thermal 
conductivity  data. 

The  fraction  of  water  removed  in  the  vapor  phase,  §,  equals  the 
ratio  of  the  vapor  flux  to  the  total  water  generation  rate,  which 
is  i/2F.  The  maximum  flux  that  can  be  removed  from  the  cathode 
through  the  vapor  phase  (calculated  in  Eq.  (23))  in  a  GDL  is  higher 
than  the  water  production  rate  if 

£  =  «(Q£  -  <2&_>v  “if)  >  1  (27) 

where  (Ls  +Lm)/D  =  (Ls/Ds)  +  (Lm/Dm),  (Ls +£m)A  =  (W^m)M£sAs)- 
D  and  A  are  the  effective  vapor  diffusion  coefficient  and  thermal 
conductivity. 

The  fraction  of  water  removed  through  the  vapor  phase  (calcu¬ 
lated  in  Eq.  (27))  determines  the  water  transport  scenario.  Different 
scenarios  for  Type  1  and  2  GDLs  are  defined  in  Table  1.  Table  2 
summarizes  the  criteria  determining  the  different  scenarios.  In  the 
water  transport  scenarios  (see  Table  1),  1  corresponds  to  under¬ 
saturated  water  vapor  at  the  cathode/MiPL  interface;  otherwise,  2 
means  saturated  vapor  at  this  interface.  In  the  parenthesis,  a  and 
b  correspond  to  under-saturated  and  saturated  water  vapor  at  the 
substrate/MiPL  interface,  respectively.  Therefore,  the  name  of  water 
transport  scenario  provides  vapor  saturation  details  in  a  GDL. 

For  Type  1  GDL,  we  predict  Scenario  1(a)  if  criterion  (27)  is  valid 
and  Scenario  2(a)  otherwise.  In  Scenario  1(a)  (Fig.  3a),  water  vapor 
in  both  substrate  and  MiPL  is  under-saturated;  therefore  100%  of 
the  water  is  removed  in  vapor  phase.  In  Scenario  2(a)  (as  shown 
in  Fig.  3b),  liquid  water  streams  in  the  GDL  are  likely  due  to  possi¬ 
ble  liquid  water  transport  from  the  MiPL/cathode  interface,  where 
water  vapor  is  in  equilibrium  with  the  liquid.  However,  in  this  case 
only  water  evaporation  occurs  in  the  GDL,  since  the  gas  is  under¬ 
saturated  in  the  GDL  bulk. 


Table  4 

Laser  flash  versus  modified  substrate  method  (MTM187)  results  [33]. 


Laser  flash  thermal 
conductivity  (W  (m  K)_1 ) 

Through-plane  thermal 
conductivity  (' W  ( m  K)_1 ) 

Thin,  low  Teflon  content 

0.035 

0.097 

MiPL  only 

Thin,  high  Teflon  content 

0.057 

0.097 

MiPL  only 

Untreated  Toray 

1.06 

0.798 

Scenarios  for  Type  2  GDL  are  more  complicated  (Table  1 ).  The 
case  with  lower  actual  WVPP  than  the  saturated  values  in  the  GDL 
belongs  to  Scenario  1(a),  as  indicated  in  Fig.  4a.  If  the  criterion  (27) 
is  not  valid  for  the  substrate,  but  valid  for  the  MiPL,  100%  of  the 
generated  water  is  removed  in  the  vapor  phase  through  the  MiPL, 
but  some  fraction  condenses  in  the  substrate.  This  fraction  should 
be  removed  through  the  liquid  phase,  corresponding  to  Scenario 
1  (b)  as  shown  in  Fig.  4b.  If  criterion  (27)  is  not  valid  for  the  MiPL,  less 
than  100%  of  water  generated  in  the  cathode  is  removed  in  the  gas 
phase.  The  vapor  in  the  MiPL  is  saturated,  and  vapor  condenses  in 
the  GDL.  This  corresponds  to  the  Scenario  2(b)  in  Fig.  4c.  Therefore, 
the  main  difference  between  Scenarios  1(b)  and  2(b)  is  the  water 
saturation  in  the  MiPL  (in  Table  1 ):  the  former  has  under-saturated 
water  vapor  in  the  MiPL  and  the  latter  does  not. 

3.  Results  and  discussion 

In  order  to  classify  GDLs  into  the  types  and  water  transport  sce¬ 
narios  that  have  been  defined,  GDL  thermal  properties  are  needed. 
GDL  thermal  conductivities  are  currently  available  through  both 
experimental  measurements  and  estimates  in  the  literature  [  18,31  ]. 
Tables  3  and  4  show  the  thermal  conductivities  of  substrate,  MiPL 
and  different  composite  combinations,  determined  by  laser  flash 
method  [32]  and  a  modified  substrate  hot  plate  method  [33].  In 
the  laser  flash  method,  the  measurement  of  the  thermal  diffusiv- 
ity  is  usually  carried  out  by  rapidly  heating  one  side  of  a  sample 
using  laser  energy  and  measuring  the  temperature  rise  on  the 


Table  3 

Laser  flash  thermal  conductivity  measurement  results  [32]. 


Sample 

Thickness,  I  @  25  °C 
(mm) 

Bulk  density,  p@  25  °C 

(gem-3) 

Temperature  (°C) 

Specific  heat, 
CpO(gK)-1) 

Conductivity,  A 
(W(mK)-1) 

Teflon-treated  Toray  paper 

0.19 

0.62 

65 

0.922 

2.13 

Thick,  high  Teflon  content  MiPL  +  Teflon-treated  Toray 

0.25 

0.61 

65 

0.870 

0.310 

Thin,  high  Teflon  content  MiPL  only 

0.045 

0.59 

65 

1.07 

0.057 

Untreated  Toray  paper 

0.20 

0.41 

65 

0.756 

1.06 

Thin,  low  Teflon  content  MiPL  +  untreated  Toray 

0.25 

0.50 

65 

0.862 

0.214 

Thin,  low  Teflon  content  MiPL  only 

0.070 

0.47 

65 

0.887 

0.035 
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Table  5 

Scenarios  of  water  transport  for  different  sets  of  GDL  parameters  @  1 A  cm-2  (indexes:  S,  substrate;  M,  MiPL;  a,  heat  fraction  removed  through  cathode  WTP). 


MiPL 

a 

Parameters 

Criteria  (17) 

Criteria  (27)  for 
the  substrate 

Criteria  (27) 
for  the  MiPL 

Scenario 

Thin,  low  Teflon  content 

0.36 

A.s  =  1.06  x  10-2  W(cmK)-1,  A.m  =  0.035  x  10-2  W(cmK)-1, 

£m  =  0.78,  IM  =  25  fxm 

False 

False 

True 

1(b) 

Thick,  low  Teflon  content 

0.22 

As  =  1.06  x  10-2  W(cmK)-1  AM  =  0.035  x  10“2  W(cmK)-1, 

£m  =  0.78,  IM  =  75  fxm 

False 

False 

True 

1(b) 

Thin,  high  Teflon  content 

0.46 

As  =  0.8  x  10-2  W  (cm  K)"1 ,  AM  =  0.097  x  10“2  W  (cm  K)-1 , 

£M  =  0.72,  IM  =  25  |Jim 

False 

False 

True 

1(b) 

Thick,  high  Teflon  content 

0.37 

As  =  0.8  x  10“2  W  (cm  K)-1 ,  AM  =  0.097  x  10“2  W  (cm  K)-1 , 

£m  =  0.72,  IM  =  75  [xm 

False 

False 

True 

1(b) 

opposite  side.  The  time  that  it  takes  for  the  heat  to  travel  through 
the  sample  and  cause  the  temperature  to  rise  on  the  rear  face 
can  be  used  to  determine  the  through-plane  thermal  conductivity. 
In  the  other  method  for  substrate  conductivity,  about  20  sam¬ 
ples  of  substrate  are  stacked  under  690  kPa  (100  psi),  and  the  heat 
flux  through  the  stack  is  measured  with  a  temperature  gradient 
applied  through  the  stack.  The  MiPL  thermal  conductivity  was  cal¬ 
culated  by  means  of  the  inverse  resistance  law  from  the  measured 
substrate  and  substrate/MiPL  composite  thermal  conductivity.  The 
measured  thermal  conductivities  of  wet-proofed  and  untreated 
Toray  substrate,  Toray/MiPL  composites  and  different  MiPLs  sur¬ 
prisingly  show  that  the  measured  MiPL  thermal  conductivity  was 
about  30%  higher  than  the  one  calculated  by  means  of  the  inverse 
resistance  law  from  the  measured  substrate  and  substrate/MiPL 
composite  thermal  conductivity.  This  may  be  due  to  the  fact  that 
the  MiPL  penetrates  into  the  substrate.  The  thermal  conductivity 
results  in  [33]  are  about  equal  for  the  untreated  Toray  substrate 
and  3  x  and  2  x  larger  than  those  in  [32 ]  for  the  low  and  high  Teflon 
content  MiPLs,  respectively.  This  may  be  attributed  to  considerable 
thermal  contact  resistance  between  the  stacked  samples.  Therefore, 
the  thermal  conductivity  data  measured  by  the  laser  flash  method 
were  used  mostly  in  this  work. 

An  analysis  of  the  experimental  results  listed  in  Table  3  and 
Table  4  implies  Type  2  GDL  and  Scenario  1  (b)  (see  Table  5)  prevailed 
at  1 A  cm-2.  That  means  that  condensation  is  not  expected  in  the 
MiPL  and  partial  condensation  is  expected  in  the  substrate  as  shown 
in  Fig.  4b.  Therefore,  due  to  sufficient  AT,  product  water  is  trans¬ 
ported  through  the  MiPL  only  in  the  vapor  phase  at  T=  65  °C.  Any 
water  vapor  that  condenses  in  the  substrate  is  efficiently  removed 
by  wicking  to  the  cathode  WTP.  This  is  due  to  its  lower  thermal  con¬ 
ductivity  of  the  MiPL  and  therefore  a  locally  high  temperature  zone 
is  maintained  near  the  cathode.  For  an  untreated  Toray  substrate 
thickness  of  0.2  mm  with  a  thin  and  low  Teflon  content  MiPL,  the 
cathode  temperature  is  higher  than  the  WTP  temperature  by  2.6  °C 
at  1 A  cm-2  using  the  data  in  Table  3. 

The  approach  provided  in  this  work  also  enables  tuning  water 
transport  scenarios  by  modifying  GDL  properties.  The  water  trans¬ 
port  scenario  depends  on  the  water  vapor  diffusion  coefficient, 
which  relies  on  MiPL  porosity  as  shown  in  Table  6.  The  fraction 


of  water  removed  through  the  MiPL  in  the  gas  phase  decreases 
from  100%  to  50%  with  porosity  decreasing  from  0.5  to  0.25.  If 
the  MiPL  porosity  is  0.5  then  100%  of  water  is  removed  by  vapor 
phase  (Scenario  1(b))  for  both  thin  and  thick  MiPLs.  The  results  for 
thermal  conductivity  obtained  in  Table  3  imply  that  100%  of  water 
is  removed  through  the  MiPL  in  the  gas  phase  even  for  the  thin, 
low  Teflon  content  MiPL  with  porosity  of  0.25.  The  experimental 
results  and  theoretical  models  obtained  for  the  current  UTC  WTP 
cell  design  also  indicate  that  the  water  generated  in  the  cathode 
@1 A  cm-2  can  be  totally  removed  through  the  vapor  phase.  In  this 
case,  the  water  vapor  in  the  GDL  can  be  under-saturated. 

Although  the  model  has  been  developed  exclusively  from  WTP 
cells,  its  application  is  not  limited.  The  work  done  by  Owejan  et  al. 
at  GM  [13]  are  a  suitable  experimental  verification  of  this  theory  in 
a  solid-plate  cell.  With  the  same  thermal  conductivity,  cracked  and 
crack-free  MiPLs  were  shown  to  give  the  same  performance  under 
different  humidification  levels.  This  indicates  that  water  vapor  dif¬ 
fusion  is  the  primary  mode  of  water  removal  from  the  cathode 
electrode  even  at  high  liquid  water  content  (300%  outlet  RH)  [13]. 
The  basis  here  is  the  dramatic  difference  in  thermal  properties  of 
carbon  paper  substrates  and  MiPLs.  The  thermal  conductivity  of 
MiPL  is  ca.  20  x  smaller  than  that  of  the  substrate  as  shown  in 
Table  3,  therefore  a  localized  hot  zone  is  created  near  the  electrode 
surface,  causing  an  under-saturation  environment  there. 

A  primary  application  for  this  work  is  to  derive  a  preferred 
composite  GDL  that  will  maximize  performance  at  high  current 
density.  To  sustain  enough  mass  transport,  the  major  fraction  of 
water  should  be  removed  through  the  gas  phase  and  condensa¬ 
tion  in  the  GDL  should  be  minimized.  That  means  that  the  desired 
scenario  is  Scenario  2(a)  and  desired  GDL  type  is  Type  1  if  liquid 
water  and  vapor  are  in  equilibrium  at  the  cathode-MiPL  interface. 
On  the  other  hand,  for  an  under-saturated  cathode-MiPL  interface, 
drying  out  of  the  cathode  needs  to  be  prevented  and  thereby  low 
ionomer  resistance  has  to  be  sustained.  Some  fraction  of  gener¬ 
ated  water  should  be  removed  through  the  liquid  phase  as  is  the 
case  for  scenario  1(b)  where  some  water  is  expected  to  condense 
in  the  substrate  where  the  WTP  can  evacuate  it  (liquid  water  is 
wicked  away  and  absorbed  in  the  WTP  more  readily  if  the  substrate 
is  hydrophilic).  For  high  current  density  with  sufficient  AT  across 


Table  6 

Scenarios  of  water  transport  for  different  MiPL  porosity  @  1 A  cm-2  (indexes:  S,  substrate;  M,  MiPL;  a,  heat  fraction  removed  through  cathode  WTP). 


MiPL 

a 

Parameters 

Criteria  (17) 

Criteria  (27)  for 
the  substrate 

Criteria  (27)  for 
the  MiPL 

Scenario 

Thin,  high  Teflon  content 

0.4 

As  =  0.8  x  10“2  W (cm  K)-1 ,  AM  =  0.097  x  10“2  W (cm  K)-1 , 

£M  =  0.5,  LM  =  25  fxm 

False 

False 

True 

1(b) 

Thick,  high  Teflon  content 

0.37 

As  =  0.8  x  10-2  W (cm  I<)-1 ,  AM  =  0.097  x  10“2  W (cm  K)"1 , 

£m  =  0.5,  Lm  =  75  |xm 

False 

False 

True 

1(b) 

Thin,  low  Teflon  content 

0.36 

As  =  1.06  x  10-2  W  (cm  K)-1 ,  AM  =  0.035  x  10“2  W  (cm  K)-1 , 
£M  =  0.25,LM  =  25|Jim 

False 

False 

True 

1(b) 

Thick,  low  Teflon  content 

0.22 

As  =  1.06  x  10-2  W  (cm  K)"1 ,  AM  =  0.035  x  10“2  W  (cm  K)"1 , 
£M  =  0.5,  Lm  =  75  fxm 

False 

False 

True 

1(b) 

Thick,  low  Teflon  content 

0.22 

As  =  1.06  x  10“2  W  (cm  K)"1 ,  AM  =  0.035  x  10“2  W  (cm  K)"1 , 
£m  =  0.25,  IM  =  75  [xm 

False 

False 

False 

2(b) 
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the  GDL,  all  water  for  this  scenario  would  be  transported  through 
the  MiPL  in  the  gas  phase.  At  smaller  current  density  with  lower 
AT,  a  larger  fraction  of  water  could  be  removed  through  the  liquid 
phase  in  both  the  MiPL  and  substrate.  This,  however,  does  not  neg¬ 
atively  impact  performance  because  the  total  water  flux  is  lower 
than  that  at  high  current  density.  Thus,  flooding  in  the  MiPL  would 
not  be  expected. 

4.  Conclusions 

A  mathematical  model  to  study  water  removal  modes  at  the 
cathode  in  a  porous-plate  cell  has  been  developed.  The  impact 
of  thermal  conductivity  and  diffusion  coefficients  on  water  trans¬ 
port  in  the  GDL  has  been  discussed.  Assuming  saturated  vapor  at 
both  GDL  interfaces,  it  is  concluded  that  a  composite  GDL  struc¬ 
ture  (with  MiPL)  with  low  thermal  conductivity  in  the  vicinity  of 
the  cathode-GDL  interface  helps  to  sustain  under-saturated  vapor 
in  the  GDL  bulk.  Water  transport  (the  fraction  of  water  removed  by 
vapor  phase  and  condensation  conditions  in  the  GDL  components) 
depends  on  the  combination  of  water  vapor  diffusion  coefficients 
and  thermal  conductivities  of  the  GDL  components.  Different  water 
removal  scenarios  were  described  (Table  1 )  and  the  criteria  to  pre¬ 
dict  the  water  removal  scenario  from  GDL  properties  were  provided 
(Table  2).  The  desirable  fraction  of  water  removed  by  vapor  phase 
is  about  90%  for  composite  GDLs,  corresponding  to  Type  2  and  Sce¬ 
nario  1(b)  or  2(b),  depending  on  the  MiPL  thickness  and  porosity 
(see  Table  6).  An  analysis  of  thermal  conductivity  impact  on  water 
transport  with  different  Teflon  content  in  MiPLs  predicts  that  more 
than  90%  of  product  water  (@  1 A  cm-2)  is  transported  through 
these  MiPLs  in  the  vapor  phase. 

Appendix  A 

A.l:  The  driving  force  of  the  vapor  transport  in  GDL 


Vapor 

MvO"i’P i)=nw(T) 

MT.P) 

Vapor 

flvOVP  2)=flw(^2) 

1 

1 

♦ 

Pi=P(T) 

Vapor  flux 

Ti>T2 

i 

P2=P(T2) 

T  Hw(Ti) 

l 

P,>P2 

~^2  MwO^) 

Water 

c,>c2 

Water 

Fig.  A.l.  In  both  of  the  vessels,  liquid  water  is  in  equilibrium  with  the  vapor  and 
the  vapor  partial  pressure  (concentration)  is  a  function  of  the  temperature  of  the 
vessel.  There  is  no  liquid  water  in  a  pipe  connecting  the  vessels  so  the  water  vapor 
concentration  and  the  temperature  in  the  pipe  are  independent. 


P  (or  concentration  C)  and  of  the  gas  temperature  T.  The  chemical 
potential  of  liquid  water  is  a  function  of  the  temperature  only 

Mw  =  /tw(T)  (A.3) 


In  liquid  water/water  vapor  equilibrium  system,  the  water  vapor 
chemical  potential  is  equal  to  the  liquid  water  chemical  potential. 
From  Eqs.  (A.l )  and  (A.3),  we  obtain  the  following  equation  for  the 
saturated  vapor  pressure  in  the  vessels 

Mw(7')  =  M°(r)  +  Rrin(P)  (A.4) 


Solving  it  with  respect  to  P,  we  obtain  the  following  dependence  of 
the  saturated  vapor  pressure  on  temperature 


p(T),ex 


(A.5) 


At  the  GDL/WTP  interface,  water  vapor  is  in  balance  with  liquid 
water.  The  water  vapor  chemical  potential  at  the  GDL/WTP  inter¬ 
face  is  equal  to  the  liquid  water  chemical  potential  in  the  WTP. 
However,  the  chemical  potential  of  the  liquid  water  in  the  cathode 
is  higher  than  the  chemical  potential  of  the  liquid  water  in  WTP 
because  the  temperature  of  the  cathode  is  higher  than  the  tem¬ 
perature  of  the  WTP.  This  implies  that  the  water  vapor  chemical 
potential  in  the  cathode  is  higher  than  the  water  vapor  chemical 
potential  at  the  GDL/WT  interface.  The  temperature  and  chemical 
potential  gradients  cause  the  thermo-diffusion  and  the  diffusion 
fluxes  the  GDL. 

To  understand  the  mechanism  of  the  vapor  transport  in  the 
liquid-water-free  GDL  we  use  a  simple  cartoon,  depicted  in  Fig.  A.l. 
The  cathode  and  WTP  are  shown  as  two  vessels,  partially  filled  with 
liquid  water,  the  GDL  is  shown  as  a  liquid-water-free  pipe  between 
these  two  vessels.  Liquid  water  is  in  equilibrium  with  vapor  in  both 
vessels.  The  temperature  of  the  left  hand  vessel,  T\,  is  higher  than 
the  temperature  of  the  right  hand  vessel,  T2,  i.e.,  T\  >T2.  AtP=  1  atm, 
the  water  vapor  can  be  treated  as  an  ideal  gas.  The  chemical  poten¬ 
tial  of  1  mol  of  ideal  gas  is 

HV(P,T)  =  II°(T)  +  RT  ln(P)  (A.l) 


Eq.  (A.5)  indicates  that  the  vapor  chemical  potential  and  pres¬ 
sure  in  both  vessels  depend  only  on  vessel  temperature.  Using  the 
Eq.  (A.2)  we  obtain  the  following  Equation  for  saturated  vapor  con¬ 
centration 

C(T)  -  exP^w(7~)  -  ^°(T))/Rr) 

Vapor  concentration  and  temperature  in  the  left  hand  side  vessel 
is  higher  than  that  in  the  right  hand  side  vessel.  This  results  in  the 
chemical  potential  gradient  and  in  the  diffusion  flux  as  shown  in 
following  Section. 

A.2:  Fields  law 

If  the  concentration  and  temperature  gradients  in  the  gas  are 
small,  the  diffusion  flux,  j,  of  one  of  the  gas  components,  such  as 
water  vapor,  is  a  linear  function  of  the  gradient  of  the  chemical 
potential  of  this  component  and  of  the  temperature  gradient  in  the 
gas [34] 

j  =  -a-  V/Xj  -  p  •  VT  (A.7) 

Eqs.  (A.4)  using  (A.2)  can  be  re-written  in  the  form 


where  P=P/Po  is  the  normalized  pressure,  p0  is  the  partial  pressure 
of  gas  under  the  standard  condition,  /z°(T)  is  standard  chemical 
potential  of  the  gas.  The  concentration  of  a  component  of  the  ideal 
gas  mixture  is 


C  = 


P 

RT 


(A.2) 


Here  P  is  the  partial  pressure  of  a  gas  component.  The  chemical 
potential  of  the  gas  component  is  a  function  of  its  partial  pressure 


Mw(T)  =  mS(T)  +  RT  In  (CjRT) 


(A.8) 


where  the  local  water  vapor  chemical  potential  /xz  is  a  function  of 
the  local  water  vapor  concentration  Q  and  of  the  local  gas  temper¬ 
ature  T.  The  gradient  of  the  chemical  potential  is  a  function  of  the 
concentration  and  temperature  gradients 


Vfc  = 


VQ  + 

T 


VT 

c 


(A.9) 
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Substituting  (A.9)  into  (A.7),  we  obtain  the  following  Equation  for 
the  flux 


j  =  -a 


VT 


(A.10) 


Using  the  standard  definitions  of  the  transport  rate  coefficients 


9T  J 


+  P 


we  finally  obtain  the  following  transport  rate  equation 


(A.11) 

(A.12) 


J  =  -d(vQ+^Vt)  (A.13) 

Here  kj  is  a  thermo-diffusion  ratio.  The  diffusion  coefficient  D 
depends  on  the  concentrations  of  all  component  of  the  gas  and 
on  the  temperature.  However,  if  the  concentration  and  temper¬ 
ature  gradients  are  small  (as  in  our  case)  it  can  be  treated  as  a 
constant. 


Having  in  mind,  that  the  water  vapor  is  saturated  in  the  cathode  and 
at  the  cathode/WTP  interface,  we  obtain  the  following  equation  for 
the  pressure  gradient 


Pc-Po  .9Po  T  VT 
lp0  ar  p0  2  t 


(A.20) 


Taking  advantage  of  Eq.  (A.20),  we  obtain  the  following  equation 
for  the  ratio  of  the  diffusion/thermo-diffusion  fluxes 


VQ/Q  =  dP0T 
VT/2T  dT  P 

At  65  °C, 

VQ/Q  _  2  x  104  x  340 
Vf/2 f  "  25  x  104 


(A.21) 


(A.22) 


Eqs.  (A.21)  and  (A.22)  implies  that,  under  typical  PEM  operational 
conditions,  the  second  term  (thermo-diffusion)  in  Eq.  (A.17)  is  much 
smaller  than  the  first  one  (diffusion)  and  that  the  conventional 
diffusion  equation 


j  =  DVQ 


(A.23) 


is  applicable  to  the  water  vapor  transport  in  GDL. 


A3:  Kinetic  coefficients  in  the  transport  equation 


References 


The  diffusion  and  thermo-diffusion  coefficient  cannot  be 
calculated  using  thermodynamic  methods.  Kinetic  theory  [35]  pro¬ 
vides  the  following  Equations  for  these  coefficients  in  the  gas 
phase 


v 

crt 


kj  =  QT— In 


(A.14) 

(A.15) 


where  v  is  the  thermal  velocity  of  the  gas  component  molecule,  at  is 
the  transport  collision  cross-section  of  the  gas  component  molecule 
(crt  practically  does  not  depend  on  the  velocity).  The  mean  thermal 
velocity  is  a  function  of  the  temperature  (v)  =  ( kT/mj )^2,  where 
m*  is  the  molecular  weight.  Thus,  for  the  thermo-diffusion  ratio  we 
obtain  the  following  equation 

fex  =  j  (A.16) 

Finally,  the  flux  of  the  gas  component  is 


j  =  DC, 


rvq 

l  Q 


VT' 

2T_ 


(A.17) 


A.4:  Fields  flux  versus  thermo -diffusion  in  the  GDL 


Here  we  compare  the  first  (diffusion)  and  second  (thermo¬ 
diffusion)  terms  in  the  right  hand  side  of  Eq.  (A.17)  using  the  typical 
GDL  conditions  in  PEM.  The  first  term  in  the  right  hand  side  of  Eq. 
(A.17)  is 


VCv  ^  Cc  -  C0 
Cv  ~  LC0 


(A.18) 


Here  L  is  the  GDL  thickness.  Using  the  equation  of  state  of  ideal  gas, 
Eq.  (A.2),  we  express  the  vapor  concentration  in  terms  of  the  partial 
pressure 


VCv  _  Pc-Po 
Cv  ~  LP0 


(A.19) 
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